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Processing Seismic Data 



The present invention relates to processing seismic data and, in particular, to processing 
marine seismic data that includes swell noise so as to attenuate or remove the swell 



noise. 



In a seismic survey, acoustic energy is emitted from a source and is detected by a sensor 
(generally referred to as a "receiver") located at a distance from the source. Some of the 
acoustic energy emitted by the source passes into the earth's interior and is reflected by 
geological structures within the earth. Information about the geological structure of the 
earth's interior can be derived from the reflected acoustic energy incident on the 
receiver. 



In a marine seismic survey both the acoustic source and receiver are disposed within a 
water column, such as the sea. Figure 1 is a schematic illustration of one marine 
seismic survey in which an acoustic energy source 1, such as an airgun, is suspended in 
the water column 2 from a survey vessel 3. An array of receivers 5 such as, for 
example, hydrophones, is also disposed in the water column 2, and the receivers 5 are 
suspended from a float 6 that is towed by a second survey vessel 4. Such a receiver 
array is known generally as a "towed marine" receiver array. 



When the acoustic source 1 is actuated it emits acoustic energy. Acoustic energy 
emitted by the source propagates generally downwards through the water column, 
passes into the earth's interior 7 and continues to propagate generally downwards 
through the earth's interior until it undergoes partial reflection at a geological structure 
within the earth's interior. Reflected acoustic energy propagates generally upwards 

through the earth's interior 7 and through the water column 2 and is incident on the 

I- 

array of receivers 5. Seismic data acquired by the receivers 5 following actuation of the 
acoustic source 1 are transmitted to the second survey vessel 4, for example by an 
electrical or optical link. The data may be analysed in a processor 8 located on the 
second survey vessel or may be stored for subsequent analysis. The principle of seismic 
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K^.,t thP Pftrth' s interior by suitable processing of the 
surveying is to obtain information about the earth s mienor uy 

seismic data acquired at the receivers 5. 

,„ gea^al. mo. *an one geological sm..u« «iU.. *e e^'s acU as a pa«i- 

for aoo..c ^,y. sei^lc ^ ac*d a. a «oeiv.r foUow^^^ 
Luon of *e acoustic so^ 1 *-fo. contain a number of spa„al„ cohc^ 
even, each of which conesponds «, ,cfle«ion of ^usUc cnerg, a. an assoc.a,ed 
geolo^cal feature wi*in eanh .ha. ac« as . parti- reflector of aco«ic cnew. 

one in processing towed n^rine seisn^c data is Uta. data acquired .n „ugh ^ 

conlns ntay contain swd, noise. Swell „o.se n»y he either coherent or random and 
at iow Lquencies, gener.., helow 30H.. SweU noise ge^ h, a ™^ ^ 

• 1 „5th thP desired events in the seismic data (i.e.. the events 
interferes destrucuvely with the desirea eveni^ . , ^. . o,„,,i 

arising from reflection of seismic energy by geolo^ca. features widnn the earti,,. Swd^ 
noiseLst he removed from the seismic data « some stage during the P— ^ ^ 
aata and the removal o, swell noise is one of the most important steps leadmg to flnal 
intcrptetation of .he acquired data and subsequent mappng of the earth's intenor. 

A number of prior art methods for attenuation of swell noise are known. 

J "A.itnmatic editing of noisy seismic data", 
G.A. McMechan has proposed, m AutomaUc ediung ;r 

X, I -JT ««8T5 8Q2 a9i9), a method in which the low 
Geophysical Prospecting, Vol. 37, pp875-892 (i^^^), ^ 

f vc.r,ir trace is evaluated. A trace is discarded if its low 
freauencY energy content of a seismic trace is evaiudwu. 

frequency gy ^ ^,hold. This method has the severe 

ftequency energy content exceeas a given 

disadvantage tha. the desired signal in the trace is also lost in the process. 

HTnmtan Holcomh. and Robert S. Wojslaw have proposed, in "Spatially weighted 
TstacMng, a technique for ptestaC suppression", S^iety of Explo^on 
Oeophysicsts, ppi.57-n60 (1992). a method that involves ^e 
JId replacing the affected portion of ti« trace with an in.«polation whrch makes 
^ Of the neighbouting traces. Since ti.s method requires an 
neighbouring t«ces reflector dip tnformatio. is therefore required, so that th,s metitod ,s 
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a costly and often inaccurate process. Furthermore the dip search can only provide the 
dominant dip. 

Guillaume Cambois and Jacques Frelet have proposed, in "Can we surgically remove 
swell noise ?", Society of Exploration Geophysicists. 1995. a noise fi.e prediction filter 
that is designed and applied in the f-x (frequency-space) domain. This filter is quite 
efficient at filtering out swell noise when it is random, but coherent swell noise is not 
eliminated by this method but rather is enhanced. TOs "one step ahead" predicUon 
filter is a coherent enhancement process and hence cannot separate coherent swell noise 
from a coherent signal 

Ali Ozbek has proposed, in a patent application entiUed "Adaptive Seismic Noise and 
Interference Attenuation Method" an adaptive beamforming approach for attenuating 
the various types of coherent noise encountered in seismic data acquisition. This method 
applies an adaptive (t-x) time space filter in those parts of the frequency-wavenumber 
space that contain the noise. Seismic signals are passed while coherent noise and 
interference components that are temporally and spatially not stationary are adapuvely 
filtered. This method can be quite effective but can be computationally intensive. 

LCanales has also proposed a method of attenuating swell noise, in the SGE 54* 
Annual International Meeting (1984). This is a "one step ahead" wiener prediction 
filtering process carried out in the f-x (frequency space) domain. Coherent signal and 
coherent swell noise are both enhanced. 

S SpiB has p«.posed, in "ActapUve Dip Rl^nng i. fte M Domain". SEO 63^ A^ua. 
ta^mational MeeUng (1993) a n«hod in which Ae p«di«ion fite co.r«pond.ng .0 
all the dips wifl.in the dala are cval»at«l throngh poljmotnial factodsation of the full 
prediction filter. All the dipping events, coherent noise include, within the input data are 
then solved for using a least squar. approach. This is a difficult task as aU the dips need 
to be accounted for. 
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Proprie^ry modules d,a. provide sWlar functa. » kno«„, -d a.esc generally 
involve comparing a Ume window widnn a pven ttace wift au^undtog »acea to check 
fo, fte presence of noiae, ua>^l, againa. a user-specified Uueahold value. Noise above 
,he specifed .lueshold is «nK,ved Such meft^ are, however, very semiove .0 fte 
Uueshold, and are likely eifter to remove .00 much from the data includrng some of the 
signal of interest or to leave some residual noise. ^ 

The present invenUon provides a meUtod of processing seismic data, the data being in 
.he ft^^uency-space domain, the method comprising the steps oft detetmimng a 
paction filter ftom the seismic data at . firs, frequency a. which swell norse ,s n« 
prese"^ and applying the prediction fil« to seismic data a, a se^d frequency at whrch 
Lell noise is p«sen,, thereby to attenuate swell noise in the seismic data a, the second 

frequency. 

The prediction filter separates fte swell noise ftom dre signal of in.e^. in the dau. and 
^..enua^s dre swell noise whUe firlly preserving d« sig^ of -enuon 
p^vides effective aUenuation of bod, cohere.. ^ random sweU noise wrthou. 
affecting the signal of interest. 

A second aspect of d» invention pK.vides a method of marine seisrmc surveying 
compri^ng actuating a somce .0 emit ac«.tic energy. ac,uiring seismic data rn dte 
ftequency domain; and processing dre dau aca»ding .0 a meUtod as defined above. 

A thrrd aspect of ti» invention provides an apparatus for processing seismic data, *e 
data being in d« frequency-space domain, the apparatus comprising .neans for 
determining a prediction filter from the s«smic data a. a firs, frequency a. whrch sweU 
no.se is not present-, and means for applying dte prediction filter to seismre daU a^^a 
second frequency at which swell noise is present, thereby to atten„a» swell norse m fte 
seismic data at the second frequency. 



The apparatus may contain a programmable data processor 
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an apparatus as defined above. 

.mviHes a storage medium comprising a program for a 
A fifth aspect of the mvenUon provides a storag 

. t. nrocessor of an apparatus as defined in above, and a sixth aspect 

data processor oi an controlUng a data processor to 

provides a storage medium contammg a program tor 

perform a method as defined above. 

Preferred features of the invention are defined in the remaining claims. 

. p.fe.ed embodiment of the invenfion wiU now be described in detaU by way of an 
niLative example, with reference to the accompanying figures m which. 

Fi^re 1 is a schematic illustration of a towed marine seismic survey; 

, A.i. trace in the fiequency domain, obtained by the 

Figure 2 illustrates a typical seismic data Uace, m the treq 

seismic survey of figure 1; 

3 U a no« or a n.^ ac.*g . e„.^- of ^ 

invention; 

r,^ 4(a, .0 4,0) ms^^ ^ 0' » 

5,a, .0 5(0, Hluaua. U-e .su>U of *e — a. appUed .o sin,o 
seismic data; 

6 is a sch^a^o .>o* of an appa«» «o««nS .0 .he presen. 

invention. 

f »w» r^pftivers 5 in the seismic surveying 
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, 1 12-) ^ --"^ -"-^ at «™e . = 0, so Ca. a ^acc . 

"T " wn, ^ 1 fro. BSUrea 4(.) and 50» .a., when ^ aoonsuo 

acousuc source. « acquired at the 

i. actuated, Uieie is initially, in the absence of no.se, no signal q 

soaice IS actuatea, me mce senerally arises from acousUc energy 

^iver. THe tot event in the seisnnc 

,at has travelled direct fro. the ^'^^^'^Z^ ^^'^ — <" - 
without entering the eardi's intenor. The t,« deUy 

- n:r:::i;rr rin^Li distance 

speed of propagation of acoustic «,ergy ' « ,„„„„^ 
between the acoustic source and the receiver. The du«t e 
evenu that arise from acoustic energy that has entered the earth s intenor 
undergone reflection at a geological feature within the earth. 

, fi.„„. 4(al to 5(c) arc shown as continuous lines. In practice, 
Ue seismic data trace of figures 4W to 5W ^.^ ^ ^^^^.^^ 

However, a r^eiver will generally ^^^^^ desampling 

''=^-^'-<"^*^-^n';llAl«aced»ref„recoi.stsofase„es 
interval 1/fs may be. for example. 1, 2 or 4ms. a 

rirete values for the paramet., one obtained in each sampling operation. 

^ seismic data trace in the time domain may be Fourier tr^tsf^med to prod^ a 
r : .e amplitude of .e seisn^c parameter as -J::;: 
, schematic illustration of a typical data .a. ■„ «c ^ 
obt^ned by transforming one of the data traces of ^„ ^ i„ .be 

^ domain by. for example, ^^^^ „ ^ 

r::t :z a — ™ Of f where 

transformed data trace of figure ^^^^ 
f k related to the sampling frequency fs by the Myquisi 

fnyq IS telatea to p^.,, o for convenience. As noted above, however, a 

shown as a continuous trace m Bgure 2 for ^^^^^^ 
seismic data trace in the time domain in pracuce consists of a series 
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rather than a continuous trace. The data trace in the frequency-space domain of figure 2 
will therefore also consist of a series of values at discrete frequencies fj, where the index 
j extends to the highest value J at which fj does not exceed the Nyquist frequency. 

Swell noise is a low-frequency noise, typically having a maximum frequency of around 
30 Hz. The frequency-space domain data trace of figure 2 therefore contains two 
regions - a low frequency region 9 that contains both the signal of interest and swell 
noise, and a higher frequency region 10 that contains only the signal of interest. In 
figure 2 the frequency region 9 that contains both swell noise and the signal of interest 
extends from frequency ft to frequency fw, and that the frequency band that contains 
only the signal of interest is in the range from frequency fn+i to fnyq. 

The lower frequency limit, f^ of the band containing swell noise will generally be zero. 
If fL is not zero, the data at frequencies below ft represents a combination of signal and 
noise other than swell noise. 

In the method of the invention, a signal prediction filter is determined from the seismic 
data trace of figure 2 at a first frequency at which swell noise is not present - that is, at a 
frequency of fn+i or greater. The prediction filter is then apphed to seismic data at a 
second frequency at which swell noise is present - that is, at a frequency between ft and 
fn - so as to attenuate swell noise in the seismic data at the second frequency. The 
method is preferably repeated to obtain prediction filters for one or more other 
frequencies at which swell noise is present in the seismic data, for example by deriving 
a filter from data at a third frequency (at which swell noise is not present) and using this 
prediction filter to attenuate swell noise in the seismic data at a fourth frequency (at 
which swell noise is present). In a particularly preferred embodiment, an associated 
prediction filter is determined for each frequency in the swell noise range from ft to fa, 
to enable swell noise to be attenuated from every frequency value in the swell noise 
range. 

In the ideal case, the prediction filter that was to be applied at a frequency fj would be 
derived from seismic data for that frequency. In the present case, however, this ideal 
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criterion would require that the prediction filter to be applied to seismic data at 
frequency in the range from ft to fw would be derived from the seismic data for that 
frequency - but this would mean deriving the prediction filter from seismic data that 
contains swell noise in addition to the signal of interest. According to the invention, 
therefore, the prediction filter for a frequency component fj is derived from the seismic 

data for a frequency that is close to fj but that is outside the frequency range within 

i 

which swell noise occurs. This is a reasonable approximation, provided that the 
characteristics of the signal of interest are identical, or very similar, for the frequency fj 
and for the frequency at which the filter is derived. 

In principle, a single prediction filter determined from data for a single frequency at 
which swell noise is not present could be applied at every frequency at which swell 
noise is present. This however has the potential disadvantage that if the filter is 
conupted for any reason, then the output obtained by applying the filter will be noisy at 
every frequency in the swell noise range. It is therefore preferable that a plurality of 
prediction filters, each determined from data at a different frequency, are used so that 
the effect of one corrupted filter will be reduced. It is particularly preferable that each 
prediction filter is applied to data at only one frequency in the swell noise range. This 
means that the number of different prediction filters required is equal to the number of 
frequencies at which swell noise is to be attenuated - and, as noted above, each filter is 
preferably determined at a frequency close to the swell noise range. 

In a particularly preferred embodiment of the invention, the filter to be applied to the 
data at frequency ft is derived from data at frequency fn+i. The frequency fn+i is, as 
noted above, outside the range at which swell noise occurs, and is the lowest frequency 
outside the range at which swell noise occurs. The signal characteristics at frequency 
fn+i should be similar to those for frequency ft. 

In this embodiment, the filter to be applied at frequency ft+i is derived from the seismic 
data at fn+i- Again, the fi-equency fH+2 is outside the frequency band within which swell 
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of interest at ft+i- 

F„rt,en«ore, in this embodimen., ,or c.h M»ency U wi*in *e f„ ^ for 
1, SV.CU noi. occur. *e a.socia«d p^dicdon fl.»r .s de»n.„ed fron, *c « 

r.. a « " " ™ 

f„ in tf« swell noise band is derived from d,e seismic im at ftequency faH-u- 

The frequency range in which .he prediction m»r .s derived *us ^ a, f„.,. and *is 
^ Jcy fulfils »o important cri.^. Fi.dy, data ac,u,red at this ftequency. or a. a 
M^er ft^nency, is assumed .0 be ftee of swell no.se. Secondly, this frequency ^ . 
Z as^ssible to frequency f. so that the signal o, intetest should be smnlar M 
UenciL f.. and f. Furihcnnore. in -Ms embodiment each filter is apphed only to 
aata at one frequency, so that the effect of a corritpted filter is nnnimtsetL 

Figure 3 is a block flow diagram of one embodiment of the pt^en, invention. 

initially a. step 1 seismic daU are obtained, for example, using a «.wed marine seisntic 
imuaiiy « » r . , . . Tk» srismic data wil consist of a 

• surveying arrangement of the type shown m figure 1. Ue setsmtc <teta 

number of uaces genetaUy similar to those shown in figure 4(a) or 5(a). 

the invention may alternatively be applied U. pre-existing seisnuc data. In this case, 
step 1 is replaced by step 2 of tetiieving pre^isting seismic dau from storage. 

A. 3 ti,e seismic data a« Fourier transformed to produce seismic d«a traces of the 
generl t;pe shown in figure 2 ^ contain information as to tite ft»^ency depended, 
of the amplitude of a seisnuc parameter. 

U should be noted titat, tn an embodiment whe,. stored seismic data ^ retrieved f^om 
storage in step 2, step 3 may be omit^d if tite retrieved seismic data are abeady « *e 
ftequency domain; Utis is indicated by a broken line. 
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At Step 4, a frequency counter / is initialised, and at step 5 a filter Fs is determined for 
an initial frequency component f,- in the range where swell noise is present in the seismic 
data. Conveniently, step 5 is initially performed for the lowest frequency component at 
which swell noise is present, that is, for frequency f, = ft (and in this case step 4 will 
have set the counter i equal to L). As mentioned above, in one preferred embodiment 

the prediction filter for frequency ft is determined from seismic data at frequency fn+i- 

! 

At step 6, the data at frequency f, for all traces, are operated on by the filter Fs 
determined for frequency f, in step 5, according to equation L 

At step 7, it is determined whether the filter has been applied a pre-determined number 
(n) times. The pre-determined number n may be specified by the user for processing a 
particular data set. If step 7 yields a "no" determination, at step 8 the filter Fs is re- 
applied to the results obtained in step 6. Step 7 is then repeated, and if a "yes" 
determination is still not obtained, steps 8 and 7 are repeated until a "y^s" determination 
is obtained. 

For example, a user could specify that the filter is to be applied three times (n = 3), in 
which case step 8 would be carried out twice (so that the filter was applied once to the 
initial data and twice to the output of the previous iterations.) If desired, the method 
may comprise the step (not shown) of a user defining the number of times the filter 
should be applied. 

In an alternative embodiment, step 7 is replaced by a quality control step (not shovwi), 
so that the filter is applied until a desired result is obtained, and is not applied a pre- 
determined number of times. For example, the residual obtained in one step may be 
compared with the previous residual. 

At step 9, it is determined whether the counter indicates that the process has been 
carried out for all frequencies in the range where swell noise is present. If the value of 
the counter i is not equal to H so f, is not equal to fn, step 9 yields a "no" determination, 
the frequency counter / is incremented at step 10, and steps 5-8 are then repeated for the 
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1. f Assuming that a "no" detennination is still 
next freauency - in this example, fui- Assuming uiai 

iauC »nm a "ye." d.— „ . oUai„«i in s.p . A "yes" —on w. b. 
„l„ed When a n,.*»i has been c«ded oa. in all frequencies in range f. .0 f„ - 
Aat is, te all frequencies at «Wch the data contains sweU no.se. 

„el,. the n^thod ntay comprises perfonning steps 5 to 8 inittal^, for f. in 
*ch case step 10 composes dec.»enang the counter. In thrs vanaUon a yes 
^nation is obt^ned at step 9 once the frequency has been reduced to f. 

O^e a "yes" deter^nation is obtained a. step 9, the results of dre n»thod .^y be 
!2 if desired (no, shown,. Tire results output ntay include, for each on^nal dau 
Z\ corrected version of ^ data trace fron, which swell noise has been re.o ed^ 
ZZ^ d.a .ace may be obtained, for exan^le, by subtracUng Ute fin^^^s^ 
trained at d« end of the iterafon st^ pocess 7 a.^ 8 f.»n. the ongtnal data t^ 
tllly or aitemauvely, the ..ults output nray include the .sidt. obta-ned f.n> 



each data trace. 



Aadidonally or al.emaa.ely .1. corrected data traces .ay be subjected » 

processing steps, to obtain infom.tio„ about the geolo^cal ^<^^^^^ 

J . * 1 o Tf Tipressarv the corrected data traces are 
inierior and this is indicated at step 12. If necessary, me i. 

intenor, ana i .,.„„„ hefoie step 12 is applied. These further 

iransfonned back to the Ume domain at step 11 before step pp 

.lata processing steps are conventional, and will not be described here further. 

Hgures 4(a, to 4(c) and 5(a) to 5(c) i,lus.«te results of the present invendon. Fig«.s 
4(c) illustrate dre invention as apphed to syntheHc sei^r^c data, and figures 5(a) 
,0 5(c) illustrate the invenUon as apphed to actual marine seismrc data. 

4(a) shows syntheUc marine seisnnc data obtained using a simuladon method 
2 ,1=1 swell noise. V. data were Simula^ for a seismic su^eytng 
tl^Lnt in Which the direct anival occurs approximately O.B seconds aft. 
rlsouree has been actuate. Frgure 4(a, shows approximately 75 simulated data 
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. •„ Each trace Shows the amplitude (horizontal axis) Of the 

.aces. . the ^ ^^^^^^^ the actuation of the acousUc 

pressure as a function of the time ( ^^^^^^ ^^^.^^^^ 

source. It will be noted that the data traces contain signitic 
times before the direct arrival. 



. figure 4(a) after sweU Doise has been attenuate using a 
„«U,od of U« invention, and figure 4(c) sho» _^ 

r::::":::::: . - ..a... .e o... ^ 

traces.) 

i fh.re is very little structure in the seismic 
- A ;f M/iil he seen that there is vciy 

"^Tr""Ie ail lvl «c« a. a^«-e., 0.S .ccona. after Ute 
aa.a tt«=es before *e d.rec< am ^ ^ ,„ 

^„^g s»eH noi. Since, in .e ab^o^-*^ 

aco^ticenergyfion-O-esourcebefo^ - ^ ^ ^ 

U significantly less sttucune ,n *e tWs indicates W 

' 3 seeds ^ - Itr: ,.re .a) a. 3 seconds or 

„„sl of d.e strucmre seen m *e onginal ^^^^^ ^ 

^,er after 0» acm«ion of *e acousdc source does no. .esul 
geological structure but arises fKm sweU noise. 

in the seismic data traces occurnrg at tijnes from 
I, wiU furdter be noted that the event, m the setsnuc ^ 

, 1 ,i 1 seconds after the actuaoon of the seismic so 
approximately 1.6 to 3 seconds ,h,„ i„ the original data traces 

dearer in the corrected seismic data traces of fi^ m 

. figure 4(.). TMs in.cates ^ of the 

^oughout the enure data trace As a r«^.^^ ^ ^^^^ ^ „ ^ 

earth'sintedorderivedfromanalysisoftheseeventssn 

. „ ^ff;mirp4ra^ had been used, 

original uncorrected data traces of figure 4(a^ 
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Figures 5(a) to 5(c) again show results of the present invention, but applied to actual 
seismic data. These seismic data were acquired using a marine seismic surveying 
arrangement in which the direct arrival occurs approximately 1.8 seconds after the 
actuation of the acoustic source. The figures show data for approximately 70 receivers, 
and these are arranged as a linear array extending generally away from the source. As a 
result, the arrival time of the direct event increases gradually from one trace to the next, 
as a consequence of the increasing distance between the source and the receiver used to 
acquire a particular trace. 

As is clearly seen in figure 5(a), the acquired seismic data contain significant amplitude 
variations for times before the arrival of the direct event and, as explained above, these 
amplitude variations arise primarily from swell noise. 

Figure 5(b) shows the traces of figures 5(a) after attenuation of swell noise using a 
method of the present invention. It will be seen that the amplitude variations in the 
traces for times before the arrival of the direct event have been significantly reduced. It 
will also be seen that many of the events in the traces for long arrival times have also 
been eliminated, and this indicates that these events arise from swell noise rather from 
reflection by a geological feature within the earth. The events of interest in the data 
traces, which generally occur at times between 1.8 and approximately 5-6 seconds after 
the actuation of the source are much clearer in figure 5(b) than in the original traces of 
5(a). 

Figure 5(c) again shows the residual, which represents swell noise. The results for the 
corrected data traces and the swell noise have again been transformed back to the time 
domain to enable easy comparison with the original data traces. 

A preferred embodiment of the filtering method will now be described. 

Let denote the input data. This is assumed to consist of N traces identified by the 

index / where 1 < / < N. These traces are acquired at the same time but by different 
receivers in a receiver array. Each trace will contain a signal of interest, and will also 
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contain swell noise. The signal of interest in the trace will be denoted by X^- and the 
swell noise value in the ^ trace will be denoted by X^i . Thus: 



The signal prediction filter is denoted by Fs = [pg, Fjj . • . Fjl], where ¥sk 

denotes the filter component. Fs is a general term for a prediction filter derived at a 
particular frequency outside the swell noise range. The index k satisfies 0 < k <P, 
where P is the number of signal related dips. The filter length is preferably between 3 
and 10 traces. The prediction filter may be derived by any suitable technique such as, 
for example, Wiener filtering. As noted above, in a preferred embodiment the 
prediction filter for ft^equency is determined from data at frequency fn+i, the 
prediction filter for frequency fL+i is determined from data at frequency fn+a, and so on. 

An appropriate filter Fs is applied to each seismic data trace at each frequency in the 
range (ft - fn) for all traces. The resultant signal is predictable as a function of space 
(Canales 1984, above). 

Filtering the input data with F^ only will result in the following system of linear 
equations: 
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The right hand side of equation (1) contains two tenns. The first term represents the 
predicted signal for each data trace. Equation (1) assumes that the filter ¥$ is exactly 
correct, so that this term is the exact signal component Xsi of each data trace D,. 

The second term on the right hand side of Equation (1) is the results of the convolution 
of the swell noise with the signal filter, and this is termed the "residual''. If the filter Fs 
is exactly correct, the residual will not contain any of the signal of interest. 

In practice, the filter Fs derived from the data traces D, is unlikely to be exact so that the* 
residual of the right hand side of equation (1) after a single application of the filter will 
contain some components of the signal of interest. It is therefore preferred to adopt an 
iterative approach, so that the residual term can be attenuated through further iterations 
thereby reducing the signal component of the residual. This is done by re-applying the 
filter Fs to the full right-hand side of equation (1), and adding the results to the 
corresponding results obtained by the first iteration. The summation of the results of the 
two iterations will enhance the quality of the signal (the terms in Xs) and attenuate the 
noise (the term in Xn) through destructive interference. Further iterations may be 
applied as necessary, until the noise has been attenuated to a satisfactory level 

Equation (1) is performed at each frequency in the swell noise range, using the 
appropriate Fs (i.e. Fs(fL) for data at frequency fu Fs(fL+i) for data at frequency ft+u 
etc). 
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If the lower frequency limit of the swell noise range. Ft, is greater than zero, data at 
frequencies below Fl will not be filtered. In this case, the output will consist of: 

a) for frequencies of Fh+i or above - unchanged signal-only data (since there is no 
swell noise in this range in the initial data); 

b) for frequencies in range from Fl to Fh - signal -only data (obtained by removal of 
swell noise); and ' 

c) for frequencies below Fl - signal and other noise (not swell noise). 

If the lower frequency limit of the swell noise range, Fl, is zero, the output will consist 
only of (a) and (b) above. 

Figure 6 is a schematic block diagram of an apparatus 11 that is able to perform a 
method according to the present invention. The apparatus 1 1 may be used in a marine 
seismic surveying arrangement in place of the conventional processor 8 of Figure 1, or 
it may be used to process pre-existing data. (If the apparatus 11 is used in a 
conventional marine seismic surveying arrangement of the type generally shown in 
Figure 1, the seismic surveying arrangement requires means for acquiring seismic data 
in the frequency-space domain and these are conveniently provided by receivers for 
acquiring seismic data in the time domain and transform means for transforming the 
data to the frequency-space domain. The transform means may be incorporated in the 
apparatus 1 1 or they may be additional to it.) 

The apparatus 11 comprises a programmable data processor 12 with a program memory 
13, for instance in the form of a read only memory (ROM), storing a program for 
controlling the data processor 12 to process seismic data by a method of the invention. 
The apparatus further comprises non-volatile read/write memory 14 for storing, for 
example, any data which must be retained in the absence of a power supply. A 
"working" or "scratch pad" memory for the data processor is provided by a random 
access memory RAM 15. An input device 16 is provided, for instance for receiving 
user commands and data. One or more output devices 17 are provided, for instance, for 
displaying information relating to the progress and result of the processing. The output 
device(s) may be, for example, a printer, a visual display unit, or an output memory. 
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Sets of seismic data for processing may be supplied via the input device 16 or may 
optionally be provided by a machine-readable data store 18. 

The results of the processing may be output via the output device 17 or may be stored. 

The program for operating the system and for performing the method described 
hereinbefore is stored in the program memory 13, which may be embodied as a 
semiconductor memory, for instance of the well known ROM type. However, the 
program may well be stored in any other suitable storage medium, such as a magnetic 
data carrier 13a (such as a "floppy disk") or a CD-ROM 13b. 



